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Abstract 

An analysis of the sensitivity of supersymmetric dark matter to variations in the b 
quark mass is given. Specifically we study the effects on the neutralino relic abundance 
from supersymmetric loop corrections to the mass of the b quark. It is known that these 
loop corrections can become significant for large tan f3. The analysis is carried out in the 
framework of mSUGRA and we focus on the region where the relic density constraints 
are satisfied by resonant annihilation through the s-channel Higgs poles. We extend the 
analysis to include CP phases taking into account the mixing of the CP-even and CP- 
odd Higgs boson states which play an important role in determining the relic density. 
Implications of the analysis for the neutralino relic density consistent with the recent 
WMAP relic density constraints are discussed. 
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1 Introduction 



There is now a convincing body of evidence that the universe has a considerable amount 
of non baryonic dark matter and the recent Wilkinson Microwave Anisotropy Probe 
(WMAP) data allows a determination of cold dark matter (CDM) to lie in the ranged |2] 

n CDM h 2 = o.ii2Qt° Z (i) 

One expects the Milky Way to have a similar density of cold dark matter and thus there are 
several ongoing experiments as well as experiments that are planned for the future for its 
detection in the laboratory [31 IH 03 EH E|- One prime CDM candidate that appears naturally 
in the framework of SUGRA modelsjS] is the neutralino We will work within the 
framework of SUGRA models which have a constrained parameter space. Thus without 
CP phases the mSUGRA parameter space is given by the parameters mo, mi , A , tan/3 
and sign(fi) where mo is the universal scalar mass, mi is the universal gaugino mass, Aq 
is the universal trilinear coupling (all given at the grand unification scale Mq), tan/3 = 
(H 2 ) I {Hi) where H 2 gives mass to the up quark and Hi gives mass to the down quark and 
the lepton, and \x is the Higgs mixing parameter which appears in the super potential in 
the form \iHiH 2 . SUGRA models allow for nonuniversalities and with nonuniversalities 
the parameter space of the model is enlarged. Thus, for example, SUGRA models with 
gauge kinetic energy functions that are not singlets of the 577(3) x SU(2) x U(l) gauge 
groups allow for nonuniversal gaugino masses (i=l,2,3) at the grand unification scale. 
The parameter space of SUGRA models is further enlarged when one allows for the CP 
phases. Thus in general /x, A and m 8 become complex allowing for phases M , a^ , and 
£j where 8^ is the /i phase, «a is the A Q phase and & is the phase of the gaugino mass 
fhi (2=1,2,, 3). Not all the phases are independent after one performs field redefinitions, 
and only specific combinations of them appear in physical processes |lUj. 

In most of the mSUGRA parameter space the neutralino relic density is too large. 
However, there are four distinguishable regions where a neutralino relic density compati- 
ble with the WMAP constraints can be found. These regions are discussed below. (I) The 
bulk region: This region corresponds to relatively small values of m and mi and is dom- 
inated by sfermion exchange diagrams. However, it is almost ruled out by the laboratory 
experiments. (II) The Hyperbolic Branch or Focus Point region (HB/FP)jTTj: This region 
occurs for very high values of rriQ and small values of /x and is thus close to the domain 
where the electroweak symmetry breaking does not occur. Here the lightest neutralino 
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has a large Higgsino component, thereby enhancing the annihilation cross- section to gauge 
boson channels. Furthermore, chargino coannihilation contributes as the chargino and the 
lightest neutralino are almost degenerate. (Ill) The stau coannihilation region: In this 
region rrif ~ m x and the annihilation cross-section increases due to coannihilations \f\. 
(IV) The resonant region: This is a rather broad region where the relic density constraints 
are satisfied by annihilation through resonant s-channel Higgs exchange. In this work we 
will mainly focus on, this resonant region. 

While there are many analyses of the neutralino relic density there are no in depth 
analyses of its sensitivity to the b quark mass. One of the purposes of this analysis is 
to investigate this sensitivity. Such an analysis is relevant since experimentally the mass 
of the b quark has an error corridor, and secondly because in supersymmetric theories 
loop corrections to the b quark mass especially for large tan/3 can be large and model 
dependent |12j. Recently, a full analysis of one loop contribution to the bottom quark mass 
(nib) including phases was given |T3] and indeed corrections to m& are found to be as much 
as 50% or more in some regions of the parameter space. Further, mj corrections are found 
to affect considerably low energy phenomenology where the b quark enters [TH II 5j. As 
noted above the m& corrections are naturally large for large tan j3 which is an interesting 
region because of the possibility of Yukawa unification^! E| arid also because it leads to 
large neutralino-proton cross-sections [T£] which makes the observation of supersymmetric 
dark matter more accessible. However, we do not address the issue of Yukawa unification 
or of neutralino-proton cross-sections in this paper. We will also discuss the dependence 
of the relic density on phases. It has been realized for some time that large phases can 
be accommodated without violating the electric dipole moment (EDM) constraints O 
|2*UI |2~T1 by a variety of ways which include mass suppression [23]; the cancellation 
mechanism |24[ l2*K]. phases only in the third generation [2E], and other mechanisms |27j . 
One of the important consequences of such phases is that the Higgs mass eigenstates are 
no longer eigenstates of CP[2H I2H1 HH EH]- It was pointed out some time ago that CP 
phases would affect dark matter significantly in regions where the neutralino annihilation 
was dominated by the resonant Higgs annihilation |29^ IH2~] . We discuss this issue in greater 
detail in this paper. 

Since the focus of this paper is on the effects of loop corrections to the b quark mass, 
we briefly discuss these corrections. For the b quark the running mass rrib(Q) and the 
physical mass, or the pole mass Mb, are related by inclusion of QCD corrections and at 
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Figure 1: Feynman diagrams responsible for the main contribution to the neutralino 
annihilation cross-section in the region of the parameter space investigated in this analysis. 
Fig. (a) gives the s-channel Higgs exchange contribution and Fig.(b) gives the t- and u- 
channel sfermion exchange contribution to the neutralino annihilation cross-section. The 
most important decay channels for large tan f3 are for / = b,r. 



the two loop level one has 
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m b (M b ) 



(2) 



where m b (M b ) is obtained from m b (Mz) by using the renormalization group equations 
and m b [Mz) is the running b quark mass at the scale of the Z boson mass defined by 



m b (M z ) = h b (M z )-j=cos(3{\ + Am b ) . 



(3) 



Here h b (Mz) is the Yukawa coupling and Am b is loop correction to m b . Now the coupling 
of the b quark to the Higgs at the tree level involves only the neutral component of the 
Hi Higgs boson and the couplings to the i?2 Higgs boson is absent. However, at the loop 
level one finds corrections to the H® coupling as well as an additional coupling to H®. 
Thus at the loop level the effective b quark coupling with the Higgs is given by[ 



- L bbH o = (h b + 5h b )b R b L H° + Ah b b R b L H° + H.c. 



(4) 



The correction to the b quark mass is then given directly in terms of Ah b and 5h b so that 

Am b = [i?e(^)tan/? + i?e(^)] (5) 
h b h b 



A full analysis of Am b is given in Ref. ^3] and we will use that analysis in this work. 
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Figure 2: The pseudo scalar Higgs boson mass as a function of Am& for fixed values 
of tan/3 of 40, 50 and 60 when m = rrn = A = 600 GeV. 
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2 CP-even, CP-odd Higgs Mixing and b Quark Mass 
Corrections 

As already mentioned we will focus on determining the sensitivity of the relic density to 
the b quark mass in the region with resonant s-channel Higgs dominance. This region is 
characterized roughly by the constraint 

2m x ~ m,A ■ (6) 

The satisfaction of the relic density constraints consistent with WMAP in this case de- 
pends sensitively on the difference 5 = (2m x — m^) which in turn depends sensitively 
on the mSUGRA parameter space. In this context the bottom mass corrections are very 
important as the value of m\ is strongly dependent on it, as shown in Fig.Q. On the 
other hand at least for the domain where the neutralino is Bino like one finds that the 
neutralino mass m x is rather insensitive to the bottom mass correction and is almost en- 
tirely determined by mi. The resonant s-channel region is only open at large tan/3. The 
exact allowed range of tan (3 depends severely on the value of the bottom quark mass. For 
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H < the resonant region is typically open for tan (3 in the range 35-45 and for /i > for 
tan p in the range 45-55. The large tan (3 regime is also interesting for other reasons, as 
in the presence of CP violation there can be a large mixing between the CP-even and the 
CP-odd states. Moreover, the CP phases have a strong impact on the b quark mass. In 
this section we discuss the relevant part of the analysis related to these effects. It is clear 
that if the CP phases influence the resonance condition, or equivalently the ratio m x /mA, 
they will have an impact on the relic density. This ratio is affected by phases mainly 
because is strongly dependent on the bottom mass correction Am& and through it on 
the CP phases. Furthermore, the Higgs couplings relevant for computing the annihilation 
cross-section depend on the CP phases. Thus we expect the relic density to be strongly 
dependent on the CP phases. 

We begin by considering the s-channel decay to a pair of fermions, as shown in 
Fig-©( a )- The Yukawa coupling correction enters clearly here in the vertex of the neutral 



Higgs with the fermion pair. The amplitude for x{Pi)x(P2] 
Higgs mass eigenstates, H k , k = 1, 2, 3 may be written as, 

1 



f(ki)f(k 2 ), mediated by 



M[ = v(p 2 ) [S' k - ifiSfo] u( Pl ] 



iM Hk T Hk 



u(h) Cfa + iCfo-ys v(k 2 ) (7) 



where 
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2M„ 



Re(A k ), 
+ Im(A k ). 



and the parameters A k are defined by 
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[rh 2 X% + rn.XH - 2^*X; X: }. 



(8) 
(9) 

(10) 

(11) 
(12) 



Here X is the matrix that diagonalizes the neutralino mass matrix so that X T M X X = 
diag(m Xl ,m X2 ,m Xs ,m X4 ) and m Xo is the lightest neutralino. Thus is the index among 
1, 2, 3, 4 that corresponds to the lightest neutralino (later in the analysis we will drop the 
subscript on Xo an d X w iU stand for the lightest neutralino). 
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Since the CP effects in the Higgs sector play an important role in this analysis, we 
briefly review the main aspects of this phenomena. In the presence of explicit CP vi- 
olation the two Higgs doublets of the supersymmetric standard model (MSSM) can be 
decomposed as follows 



m 







f v i + 0i + iipi 


HI j 




V HT' 






( Hf 


Hi) 


V2 


\ V 2 + 02 + ilp2 



(13) 



where 0i, 02, ^1,^2 are real quantum fields and Oh is a phase. Variations with respect to 
the fields give 
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1 ,dAV, 



v 2 dipi 



(14) 



where m%, m 2 , m 3 are the parameters that enter in the tree-level Higgs potential, i.e., Vq = 
m\\H\\ 2 + m 2 l \H2 \ 2 + (m 2 ,Hi.H 2 + H.c.) + Vd where Vd is the D-term contribution, g 2 and 
gy are the gauge coupling constants for SU{2) and U{l)y gauge groups, and AV is the 
loop correction to the Higgs potential. In the above the subscript denotes that the 
quantities are computed at the point where 0i = <p 2 — ipi — ^2 — 0. Eq. (jl4|) provides a 
determination of Oh- Computations in the above basis lead to a 4 x 4 Higgs mass matrix. 
It is useful to introduce a new basis {0i, 2 , rfiD, ^zd} where ipiDi^zD are defined by 



ipiD = sin /3ipx + cos f3ifj 2 
i>2D = - cos fiipi + sin /3ip 2 



(15) 



In the new basis the field i/j 2 d exhibits itself as the Goldstone field and decouples from 
the other three fields {0i, 02, ipm} and the Higgs mass matrix in the new basis takes on 
the form 
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Higgs 



/ M^ + m^ + A 
-(Ml + m 2 A )s p cp + A12 



-(M 2 +m 2 A )s(]C f} + A 12 
M 2 z s} + m\ 4 + A 22 
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A 23 



A13 
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{m\ + A 33 ) J 



(16) 
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where is the mass of the CP-odd Higgs boson at the tree level, Mz is the Z boson 
mass, sp(cp) = sin/3(cos/?), and Ay are the loop corrections. These loop corrections 
have been computed from the exchange of stops and sbottoms in Refs. |28| I29j. from the 
exchange of charginos in Ref.j2H] and from the exchange of neutralinos in Ref.j3U]- Thus 
the corrections Ay (i,j=l,2,3) receive contributions from stop, chargino and neutralino 
exchanges. Their relative contributions depend on the point in the parameter space one 
is in. We denote the eigenstates of the mass 2 matrix of Eq.([16[) by (k=l,2,3) and we 
define the matrix R with elements Rij as the matrix which diagonalizes the above 3x3 
Higgs mass 2 matrix so that 

RM 2 iggs R T = dtag(M 2 Hi ,M 2 H2 ,M 2 Ha ). (17) 

and thus we have 







f <h \ 




= R x 








Uld/ 



In the analysis of this paper we work in the decoupling regime of the Higgs sector, charac- 
terized by rriA Mz and large tan (3. In this regime the light Higgs boson is denoted by 
H2 and the two heavy Higgs particles are described by H\ and if 3. For the case when we 
have CP conservation and no mixing of CP even and CP odd states, we denote the heavy 
scalar Higgs boson by H (at large tan/5 it is almost equal to 0i) and the pseudo scalar 
Higgs boson by A. Returning to the general case with CP phases, in the decoupling limit 
the heavy Higgs states are almost degenerate and moreover have nearly equal widths, i.e., 

m Hl ^ m H3 , ~ r# 3 . (19) 

Furthermore, the lightest Higgs boson behaves almost like the SM Higgs particle. This 
means that there may be considerable mixing between the two heavy CP eigenstates, H 
and A, whereas the mixing with the lightest Higgs is tiny. Corrections to Yukawa coupling 
arise through the parameters C'f/f that enter in Eq.© so that 

C^ k = C^cos Xb -C[ k sm Xb , (20) 

and 

Cb,k = Cb,k sin Xb + C£ k cos Xb, (21) 
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where 



V2C[ k = Re{h b +5h b )R kl + [-Im(h b +5h b ) sin f3+Im(Ah b ) cos (3}R k3 +Re{Ah b )R k2 (22) 
and where 

>/2C£fc = -Im(h h + <fa 6 )-R*i + [-Re(h + sin/? + Re(Ah b ) cosp}R k3 - Im(Ah b )R k2 

(23) 

and the angle Xb is defined by 

Im r§h k + A^tan/3) 

(24) 

The phases enter in a variety of ways in the model. Thus the parameters Qq and Rq Q 
contain the combined effects of the phases 9^, £i and £2- Similarly, Rij contain the 
combined effects of the above three phases and in addition depend on a^, (of which the 
most important is a^J. Further, C s,p derive their phase dependence through R^ and in 
addition depend on £3 which enters via the SUSY QCD corrections Ah b and 5h b . Including 
all the contributions any of the phases may produce a strong effect on the relic density. 
Explicit analyses bear this out although the relative contribution of the different phases 
depends on the part of the parameter space one is working in. The s-channel annihilation 
cross-section for x{Pi)x{V2) f{ki)f{k 2 ) is proportional to the squared of the amplitude 
given in Eq. ([?|) and reads 



M/(M/r " {-mC+ s -Mh^ (25) 

Therefore, the imaginary couplings, S^', will yield the dominant contribution to the ther- 
mally averaged annihilation cross-section, as the real couplings, S k , are p-wave suppressed 
by the factor (1 — Arrit/s). In the case of vanishing CP-phases the pseudo scalar mediated 
channel thereby dominates over the one mediated by the heavy scalar Higgs. However, 
the contribution from H mediation cannot be neglected, as its contribution is typically 
about 10%. In the presence of non-zero phases both of the heavy Higgs acquire imaginary 
coupling and both may give a significant contribution. We may neglect the contribu- 
tion from the lightest Higgs exchange diagram since it is not resonant and moreover is 
suppressed by small couplings 3 . 

3 The region where the lightest Higgs is resonant is almost excluded by laboratory constraints. 



(C^C* + C£ fc C£) \S' k S[{\ - Amy s) + Sffi 
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tan(3 

Figure 3: An exhibition of the sensitivity of the relic density to the b quark mass as a 
function of tan (3 for the case when m = mi = A = 600 GeV for values of Amj from 
(-20%) - (+20%). 

As already mentioned the inclusion of the CP phases has two major consequences; 
it affects the SUSY correction to the bottom mass Am& and it also generates a mixing 
in the heavy Higgs sector. We discuss now in greater detail the effect of mixing in the 
heavy Higgs sector. We begin by observing that in the CP conserving case the pseudo 
scalar channel gives the main contribution. As the Higgs mixing turns on the pseudo 
scalar becomes a linear combination of the two mass eigenstates Hi,H 3 , whereas H 2 stays 
almost entirely a CP-even state. However, the total annihilation cross-section which is a 
sum over all the Higgs exchanges remains almost constant. Since CP even and CP odd 
Higgs mixing involves essentially only two Higgs bosons, we may represent this mixing 
by just one 2x2 orthogonal matrix rotation. Such a rotation does not change the sum 
of the squared couplings of the two heavy Higgs boson, and thereby the effect of the 
mixing on the annihilation cross-section is small. The basic reason for the mixing effect 
being small is because of the near degeneracy of the CP even and CP odd Higgs masses 
and widths, i.e., the fact that m rriH s , and ~ r# 3 . We note in passing that 
the contribution from the Higgs exchange interference term Hh — Hf. to the neutralino 
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annihilation cross section is negligible. These phenomena allow us to write the total 
s-channel contribution in a simplified way. Thus recalling that the lightest Higgs gives 
almost a vanishing contributions, we only have to sum over the heavy Higgs particles in 
Eq . (|25|) . As we are in the decoupling limit given by Eq.(|19p. the propagators in Eq.()25|) are 
identical and can be factored out. Furthermore, for large tan/3 we have the approximate 
relations between the bottom-Higgs couplings in the CP-conserving case, 

C& * -C%, C S A ~ C£ . (26) 

These relations are independent of rotations in the Higgs sector, i.e., Higgs mixing, as is 
easily checked. Also because of the decoupling of the light Higgs boson, the mixing of the 
Higgs is described by just one angle so that 

H 1 \ = f<xx{9) sin(9)\(H\ 
H 3 J \-sm(9) cos{9)J \Aj 

which gives: 



and 



C£ = cos(0)C£ + sm(9)C s A (28) 
C^ = -sm(e)Cl+cos(6)C s A (29) 



C b P A = cob(0)C^ + sin(0)C£ = cos(9)C s A - sin^C* = Cg 3 (30) 
= - fdn(0)C£ + cos(0)C^ = - sin(0)C| - cos(0)C£ = -C£ (31) 



This is just Eq. (J26|) in the Higgs rotated basis and we see that the interference terms are 
very small 

+ C & C & * • ( 32 ) 
Furthermore, using Eq. (j25J) it is clear that the s-channel contribution is proportional to 
C s where 

C s = ((C^Y + (C£) 2 ) S? + ((C 6 5 3 ) 2 + (CQ 2 ) Sf (33) 
The b quark couplings factors out, due to Eq.(|30|). and we get, 

C s = ((Cf,) 2 + « 3 ) 2 ) (Sf + ST) • (34) 
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Again, the Higgs mixing does not change the square of the imaginary Higgs- x~X coupling. 
In the CP conserving limit we get 

5; =0, S» A = - - Qoo MP) + cot(/3) (Qgo cosOfl) + iQ (35) 

and 

(S? + 5f ) (sin(^) 2 + (cos(fl)^) 2 = (^) 2 (36) 

Thus Cs is unaffected by the Higgs mixing, but can vary with phases if the magnitudes 
and |Cf 3 | 2 + |Cf 3 | 2 vary with phases. As already discussed the CP phases have a 
large impact on the relic density through their influence on the b quark mass via the loop 
correction Am&. An exhaustive analysis of the dependence of Awij on phases is given in 
Ref.|13j. For large tan/3 and small Aq the dominant contribution to Awij comes from the 
gluino-sbottom exchange diagram and the important phases here are 9^ and £3. However, 
if Aq is large the stop-chargino correction would be large and the phase otA t plays an 
important role. There are also neutralino diagrams but normally their contributions are 
small. Thus, the CP phases 6^, £3 and aA may strongly affect the relic density, whereas 
only weak dependent on £1,^2 will be present. 

3 Sensitivity of Dark Matter to the b Quark Mass 
without CP Phases 

While a considerable body of work already exists on the analyses of supersymmetric dark 
matter (for a small sample see Ref.f|35p). no in depth study exists on the sensitivity of 
dark matter analyses to the b quark mass. In this section we analyse this sensitivity of 
the relic density to the b quark mass for the case when the phases are set to zero. In the 
analysis we use the standard techniques of evolving the parameters given in mSUGRA at 
the grand unification scale by the renormalization group evolution taking care that charge 
and color conservation is appropriately preserved (for a recent analysis of charge and color 
conservation constraints see Ref.j2E])- We describe now the result of the analysis. (For 
a partial previous analysis of this topic see Ref [37| ) • One of the parameters which enters 
sensitively in the dark matter analysis is the mass of the CP odd Higgs boson tua- Fig. (j2J 
shows rriA as a function of the b quark correction Amj, which is used as a free parameter. 
The ranges chosen are such that the tua may lie in the resonance region of the annihilation 
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Figure 4: The region allowed by the relic density constraints in the tan j3 - Am& plane for 
the case when m = mi = A = 600 GeV. Curves with fixed rriA are also shown. 
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of the two neutralinos. We find that rriA shows a very significant variation as Ami, moves in 
the range —.3 to .3. Fig. (j2J) demonstrates the huge sensitivity of to the b quark mass. 
Fig. El also shows that for fixed tan/3 one can enter in the area of the resonance for certain 
values of Ar/v Fig.|3]shows the sensitivity of the relic density to corrections to the b quark 
mass. The analysis was carried out using Micromegas^^. The dominant channels that 
contribute to the relic density depend on the mass region and are as follows: In the region 
2m x -C rriA the main channels are XX ~^ T ^ an d XX ~> bb. Here typically Q x h 2 > 0.5 
and the main contribution comes from t- and u-channel exchange of the sbottom and 
stau sparticles. Moreover, also the effects of the fi and e decay channels can be seen. 
Since their contributions are suppressed by the corresponding slepton masses, it signifies 
that one is far away from the s-channel Higgs resonances. In the region 2m x ~ the 
resonant channels account for almost the full contribution to Q x h 2 and their influence can 
be detected several widths, Fa, away from the resonance. In this region the contribution 
to the neutralino relic density from the t- and u-channel exchanges can be as much as 
10% within the relic density range allowed by the WMAP data. 

Another contribution that can potentially enter is coannihilation. Indeed for < 
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Figure 5: Region in the rriQ—m x plane of mSUGRA allowed by the relic density constraints 
for the case when tan/5 = 50, \A \ = mi for values of Amj = (left) and Ami, = 20% 
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(right). The limiting lines close areas such that Q x h 2 is below the indicated value. 

1.25xm x , one has important effects from fi% coannihilation. These effects can be observed 
at the end of the lines of Amj = 10%, 20% in Fig. (jH)). Thus for low values of tan/3 one 
is in the non resonant region and increasing tan (5 moves one to the resonant region and 
consequently the relic density decreases due to resonant annihilation. As tan (3 increases 
further, the relic density increases to become flat due to the non fermionic decays. Finally, 
the curves for Ami, of 10% and 20% exhibit coannihilation and Q x h 2 decreases again due 
to this effect as tan/? increases. In Fig. (J^) regions with fixed corridors of the relic density 
are plotted in the tan j3 — Ami, plane. The region consistent with the current range of 
relic density observed by WMAP is displayed in the dark shaded region. The hatched 
region has a value of Q x h 2 below the WMAP observation. Furthermore, curves with 
constant values of are exhibited. The analysis shows that the region consistent with 
the WMAP relic density constraint is very sensitive to the Amj correction. 

Fig.(jSJ) displays area plots in the tuq — m x plane of the relic density. The light 
shaded/grey region has 0.1291 < Q x h 2 < 0.3. The medium shaded/green region has 
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0.094 < Vt x h 2 < 0.1291 and is the 2a allowed WMAP region. Finally the dark shaded/red 
region has Q x h 2 < 0.094. The hatched region is excluded as the f is the LSP. In Fig.flHJ) 
we considered two values of Am?,; Am;, = and Ami, = 20%. A comparison of these two 
exhibits the dramatic dependence of the various regions on the b quark mass. Specifically 
it is seen that the region consistent with the WMAP constraint is drastically shifted to- 
ward lower values of m x /m for smaller Amj. Inclusion of the experimental bounds from 
processes such as b — > on Fig. (jSJ) is beyond the scope of our study. However, the case 
of Ami, = 20% is comparable to the mSUGRA case (where Amj ranges approximately 
from 17 to 21%). Therefore, the restrictions from the bounds on b — > sj, along with 
other constraints, on the left graph of Fig. JSJ) can be approximately deduced from the 
appropriate figures of Ref. [SHI- We note that in Fig.fjSJ) there is a region where m and mi 
get large and appear to have a superficial resemblance with the Hyperbolic Branch/Focus 
Point (HB/FP) region^^. However, the mechanism by which relic density constraint 
is satisfied in the WMAP region is entirely different in this case than in the HB/FP 
case (23 HOI- Thus in the analysis presented here the relic density constraints are satisfied 
by the mechanism of proximity to a resonant state (see also in this context Ref.|41|) while 
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Figure 7: A plot of {ctvm^i) as a function of £3 and 9^ (with the other phases set to zero) 
for the case when m = mi = A = 600 GeV, tan/5 = 50 and using the theoretically 
predicted value of Awij (black lines), Awij = (light lines). The contribution of dominant 
channels to (ctumcSi) are a ls° shown: all contributions (thick lines), only s-channel Hi 
mediated annihilation to bb (dashed lines) and only s-channel H 3 mediated annihilation 
to bb (dot-dashed lines) and all s-channel annihilation to rf (solid thin lines). 



for the HB/FP region the satisfaction occurs with a significant amount of coannihilation. 
In Fig.© we give a plot of Q x h 2 as a function of m x for fixed rriQ (i.e., m = 600 GeV) 
and tan/3 = 50 for Am^ values varying in the range (—10%) — (+20%). Again one finds 
that the relic density is sharply dependent on the b quark mass correction. 



4 Sensitivity of Dark Matter to the b Quark Mass 
with CP Phases 

We now give the analysis with inclusion of CP phases. In the calculation of the relic 
density, we only consider the contribution from the s-channel exchange of the three Higgs 
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yrad) yrad) 

Figure 8: A plot of the neutralino relic density Q x h 2 as a function of £3 and 9^ for 
m = mi = A = 600 GeV, tan/5 = 50 and using the theoretically predicted value of 
Arrib (black lines), Ami, = (light lines). The contribution of dominant s-channels to the 
relic density are represented by the same type of lines as fig. (J7J). 

Hi, H2, H3 and the t- and u-channel exchange of sfermions as shown in FigO The pre- 
diction for the Higgs masses and widths are extracted from the newly developed software 
package CPsuperHj32]- The impact of the CP phases on the relic density is as in the 
case without CP phases, i.e., mainly through Am&. On the other hand the effects of the 
Higgs mixing are marginal. In Fig. (J7J) we give a plot of the thermally averaged annihi- 
lation cross-section at a temperature of the order of the freeze out temperature Tf (see 
Eq.(|4U|) of the Appendix) as a function of 9 ^ and £3 for tan/3 = 50. The contribution 
of individual channels are also displayed. The channels with the bb final state dominate 
over the channels with rf final state due to the color factor. We plot Q x h 2 as a function 
of 9^ and £ 3 in Fig. (JHJ) for the same case. Figs. (JJJ) and (JBJ) also exhibit the dependence 
on the bottom mass correction, as two different values are used; the theoretical value of 
Amj, (black lines) and Awij = (light lines). The large effects of CP phases on the relic 
density in this case are clearly evident. In particular it is seen that the largest impact 
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Figure 9: Same as Fig. (JZj) and Fig. (JHJ) for tan/3 = 40. 



from the CP phases arises from their influence on the value of Am^. The curves with Amj 
confined to a constant vanishing value show less variations with the CP phases, although 
Q x h 2 still changes by almost a factor of two due to the variation of the Higgs couplings. 
The large effect of the Awij arises via its effect on vtla- Similar plots as functions of £3 
are given in Fig. (jHJ) for tan (3 = 40. 

The dependence of the neutralino relic density on aA is displayed in Fig. (|10jl. This 
dependence arises from the effect of oia on and tua- Thus for fixed A , variations 
in aA affect m fl which can generate fx coannihilations, and even push below m x . 
In Fig. (jll)) the neutralino relic density is displayed as a function of tan (3 for three cases 
given by: (i) m = m 1/2 =\A \ = 300 GeV, a Ao = 1.0, £1 = 0.5, £ 2 = 0.66, £ 3 = 0.62, 
0„ = 2.5; (ii) m = m 1/2 = |A | = 555 GeV, a Ao = 2.0, £i = 0.6, 6 = 0.65, £3 = 0.65, 
0„ = 2.5; (iii) m = m 1/2 = \A \ = 480 GeV, a Ao = 0.8, 6 = 0.4, £ 2 = 0.66, 
£3 = 0.63, 9^ = 2.5. In all cases the EDM constraints for the electron, the neutron and 
for 199 Hg are satisfied for tan/3 = 40 and their values are exhibited in tabled These 
results may be compared with the current experimental limits on the EDM of the electon, 
the neutron and on 199 Hg as follows: \d e \ < 4.23 x 10~ 27 ecm, \d n \ < 6.5 x 10 _26 ecm and 
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Case 


\d e \e.cm 


\d n \e.cm 


C Hg cm 


(i) 


2.74 x 1(T 27 


1.79 x 1(T 26 


8.72 x 1(T 27 


(ii) 


1.29 x 10~ 27 


1.82 x 1(T 27 


6.02 x 10" 28 


(iii) 


9.72 x 10~ 28 


4.19 x 10~ 26 


1.41 x 10~ 27 



Table 1: The EDMs for tan/5 = 40 for cases (i)-(iii) of text. 

C*H g < 3.0 x 10~ 26 cm from the 199 Hg analysis (where Ca g is defined as in Ref. |25j). From 
Fig. © and Fig. (j3J) it is apparent that larger negative corrections to the bottom quark 
mass push the resonance region toward smaller values of tan/3. Use of nonuniversalities 
for the gaugino masses, including the case of having relative signs among them, allows 
for larger negative corrections to the b quark mass. Therefore, it is possible to achieve 
agreement with the WMAP result for lower values of tan/3 than in the mSUGRA case. 
Considering only the main contributions from the gluino-sbottom, the bottom quark mass 
correction will reach its maximum negative value for #^+£3 = vr. The phase of the trilinear 
coupling also plays a role through the chargino loop. Thus we investigate the case with 
9^ = 0, £3 = 7r, and take oa = vr at the GUT scale. As an illustration we show in Fig. (jl2|) 
that indeed the WMAP result is compatible with tan/3 = 30 in the resonant s-channel 
region. The analysis also implies that the upper limit on the neutralino mass will be 
larger than in the mSUGRA case. For tan (3 = 30 we find an upper bound of ~ 700 GeV, 
as seen in Fig. (fT2J) . For comparison the upper bound in mSUGRA is found to be 500 
GeV for tan/? < 30 in Ref.[M]. 

5 Conclusion 

In this paper we have carried out a detailed analysis to study the sensitivity of dark 
matter to the b quark mass. This is done in two ways: by assuming that the correction 
to the b quark mass in a free parameter and also computing it from loop corrections. In 
each case it is found that the relic density is very sensitive to the mass of the b quark 
for large tan j3. In the analysis we focus on the region where the relic density constraints 
are satisfied by annihilation through resonant Higgs poles. The analysis is then extended 
to include CP phases in the soft parameters taking account of the CP-even and CP-odd 
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Figure 10: The left graph shows the dependence of (solid lines) and (dashed 
lines) on a,A a for mo = mi = 600 GeV, tan (3 = 50 and for three different values of 
|^4o|/ m i (indicated on the curves). The neutralino relic density for the same three cases 
is displayed in the graph on the right. 

Higgs mixing. Sensitivity of the relic density to variations in the b quark mass and to 
CP phases are then investigated and a great sensitivity to variations in the b quark mass 
with inclusion of phases in again observed. These results have important implications for 
predictions of dark matter in models where tan (3 is large, such as in unified models bases 
on 5*0(10), and for the observation of supersymmetric dark matter in such models. 
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APPENDIX A: Relic Density Analysis 

The analysis of neutralino relic density must be done with care since one has direct channel 
poles and one must use the accurate method on doing the thermal averaging over these 
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Figure 11: The neutralino relic density as a function of tan/3 for the three cases (i), (ii), 
(hi) of the text. Lines (I), (II) and (III) correspond to similar set of SUSY parameters for 
the case of vanishing phases. 
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Figure 12: The neutralino relic density as a function of the lightest neutralino mass for 
tan/3 = 30, £i = 0, £2 = 0, £3 = n, ctA = vr and 6* M = 0. Two different values of m are 
displayed: (i) mo = 600 GeV (solid line) and (ii) mo = 750 GeV (dashed line). 



poles [53] ■ We give here the basic formulas for the relic density analysis 



Vl x h 2 = 2.755 10 s x , 
x GeV 



x v 



The evolution equation for Y is given by 



dY 



(av Mfi )(Y 2 -Y e 2 q ), 



(37) 



(3f 



45G N "" 1V1 ~ eq> 
Here {avMfi) is the thermal average of the neutralino annihilation cross section multiplied 
by the M0ller velocity j33], Y = Y(T = T = 2.726/^), where T is the microwave 
background temperature, Y cq = Y cq (T) is the thermal equilibrium abundance given by 



F eq (T) = 2 x 



(f)V 2 (^). 



(39) 



4n*h cS (T) 

The number of degrees of freedom is g* ~ 81. However, we use a more precise value as 
a function of the temperature obtained from Ref. jlH] and the same is done for h c g. To 
calculate the freeze-out temperature Tf we use the relation 

d\n(Y eq ) 



dT 



45G 



{av M fi)Y cq 5(5 + 2) 



(40) 
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The equation for Y is 

Y = Y + XT f (41) 

and Yf = Y(Tp) = (1 + 5)Y cq (Tf). We have introduced the amount Xx f such that we can 
split the independent contribution of each channel 

Xt, = Q 9,{T)^(av m )dT (42) 

We have taken T = and 5 = 1.5 following the suggestion of Micromegas. As stated 
already care must be taken in computing thermal averaging since one must integrate over 
the direct channel poles properly [13]. We use the relation 



■ jT dsa(s)(s - AmD^sK, (43) 



8m x TK$(m x /T) 

To calculate a(s) from the partial amplitudes we use the following definitions 

a(s) = 2w(s)/^s(s-Am 2 x ) (44) 
Since we only consider channels XX //> (/ = b,r), w(s) becomes, 

= o«Z £ c / x (1 - — (45) 

OZ7T , , S 
f=o,T 

Cf is the color factor so that q, = 3, c T = 1. The definition of Wf(s) is directly related to 
the amplitude 

w f (s) = \ £ dcos0 O M|A(xx - //)| 2 (46) 
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